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Abstract 

There have been lots of efforts on the construction of quasi-cycUc (QC) low-density parity-check 
(LDPC) codes with large girth. However, most of them are focused on protographs with single edges 
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QC LDPC codes lifted from them can potentially have larger minimum Hamming distance. In this paper, 
all subgraph patterns of multiple-edge protographs, which prevent QC LDPC codes from having large 
girth by inducing inevitable cycles, are fully investigated based on graph-theoretic approach. By using 
combinatorial designs, a systematic construction method of multiple-edge protographs is proposed for 
regular QC LDPC codes with girth at least 12 and also other method is proposed for regular QC LDPC 
codes with girth at least 14. A construction algorithm of QC LDPC codes by Ufting multiple-edge 
protographs is proposed and it is shown that the resulting QC LDPC codes have larger upper boimds on 
the minimum Hamming distance than those lifted from single-edge protographs. Simulation results are 
provided to compare the performance of the proposed QC LDPC codes, the progressive edge-growth 
(PEG) LDPC codes, and the PEG QC LDPC codes. 
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I. Introduction 

Low-density parity-check (LDPC) codes [[H have been one of major research topics in coding 
area over the past decade due to their near capacity-approaching performance. Since low decoding 
complexity can be achieved by various iterative decoding algorithms, LDPC codes have been 
adopted in many practical applications. Especially, quasi-cyclic (QC) LDPC codes are well suited 
for hardware implementation using simple shift registers due to the regularity in their parity- 
check matrices. 

Thorpe [2] introduced the concept of protograph-based LDPC codes, a class of LDPC codes 
lifted from protographs. QC LDPC codes belong to the protograph-based LDPC codes because 
they can be regarded as the lifted ones from the protographs using cyclic permutations. Therefore, 
the performance of QC LDPC codes mainly depends on how to design their protographs as well 
as how to assign shift values. 

The performance of LDPC codes under message-passing algorithms depends on the girth of 
the codes because a message sent by a node along a cycle propagates back to the node itself after 
some iterations, which causes the dependence among messages and performance degradation. 
Therefore, there have been lots of efforts to construct QC LDPC codes with large girth flBll- lTm . 
In flU, necessary and sufficient conditions on determining the girth of QC LDPC codes from 
circulant permutation matrices are derived and some families of QC LDPC codes are constructed. 
Most of QC LDPC codes with large girth are constructed based on algebraic structures [|3l-[[5l, 
[l9]l- [fTT]| while some optimization algorithms and greedy search algorithms are used to find QC 
LDPC codes with large girth [|6l-[l8l. Various combinatorial designs have also been widely used 
to construct QC LDPC codes in order to guarantee the girth at least 6 [|T0l - [[T4l . 

The girth of QC LDPC codes constructed from protographs is determined by the structure of 
the protograph, the lift size, and all the shift values. In [|3|, [[TOl , and IfTSl , an upper bound 
on the girth of QC LDPC codes, which is only determined by the structure of the protograph, 
is discussed. Especially, in all substructures of multiple-edge protographs, which inevitably 
give rise to cycles of length up to 12, are searched but any construction method of multiple-edge 
protographs for QC LDPC codes with large girth is not provided. In [fTOl . all substructures of 
single-edge protographs, which inevitably give rise to cycles of length up to 20 in QC LDPC 
codes, are identified and by using combinatorial designs, some single-edge protographs for girth 



October 16, 2012 



DRAFT 



3 



larger than or equal to 18 and other single-edge protographs for girth larger than or equal to 14 
are constructed. 

Although the behavior of iterative message-passing decoders is mostly dominated by the 
pseudo-weight of pseudo-codewords [fT6l . [fTTll . the minimum Hamming distance still plays an 
important role because it characterizes the undetectable errors and provides an upper bound 
on the minimum pseudo-weight of a code. Smarandache and Vontobel fflSl derived two upper 
bounds on the minimum Hamming distance of QC LDPC codes, where one bound is applied 
when QC LDPC codes are explicitly given and the other bound can be applied even when only 
the protographs are given. It is shown by experiments that these upper bounds are very close to 
the real minimum Hamming distance when the lift size for a protograph is large enough. Also, 
through several examples, we can see that for the given protograph size and the given row- and 
column- weights, these two upper bounds increase as the number of multiple edges increases in 
the protograph. Therefore, these upper bounds can be increased if multiple-edge protographs are 
used to construct QC LDPC codes, compared to the case of single-edge protographs. 

In this paper, multiple-edge protographs which can be lifted to QC LDPC codes with large girth 
are investigated. Search for all single- and multiple-edge subgraphs which inevitably generate 
cycles of any length in QC LDPC codes are systematically performed based on graph-theoretic 
approach as an extension of the results in [fTOl . and [fTSl . Construction methods of multiple- 
edge protographs using various combinatorial designs are proposed and a lifting algorithm to 
construct regular QC LDPC codes with large girth is also proposed. 

The remainder of the paper is organized as follows. Section |ll] introduces QC LDPC codes, 
protographs, and the concept of inevitable cycles. In Section |in[ all single- and multiple-edge 
subgraphs which generate inevitable cycles in QC LDPC codes are fully searched. Based on 



these subgraph patterns. Section IV describes a design method for multiple-edge protographs 
of regular QC LDPC codes having girth larger than or equal to 12. In Section |V| construction 
methods of multiple-edge protographs are proposed for regular QC LDPC codes having girth 
14 when the variable node degree is 3 and they are generalized for regular QC LDPC codes 
with variable node degree larger than 3. In Section |VI} a construction algorithm of QC LDPC 
codes lifted from the multiple-edge protographs is proposed. It is also shown that the proposed 
QC LDPC codes have larger upper bounds on the minimum Hamming distance than those lifted 
from single-edge protographs and the performance of the proposed QC LDPC codes is verified 
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via numerical analysis. Finally, the conclusions are provided in Section |VII 



II. Inevitable Cycles of QC LDPC Codes 

A. QC LDPC Codes 

Let C be a binary LDPC code whose parity-check matrix if is a J x L array of zxz circulants 



or zero matrices as 



H 



Hi ft 



Hj-lfi -f^J-1,1 ■ ■ ■ -f^J-l,L-l_ 

where a circulant Hj^i is defined as a matrix whose each row is a cyclic shift of the row above 
it. Such an LDPC code is called quasi-cyclic because applying circular shifts to the length-z 
subblocks of a codeword gives another codeword. Also, a bipartite graph which has H as its 
incidence matrix is called the Tanner graph of C. 

The weight of a circulant Hj i is defined as the number of nonzero elements in the first column 
and denoted by wt(Hj i). A circulant is entirely described by the positions of nonzero elements 
in the first column. Let i, < i < z — 1, he the index of the (i + l)-st element in the first 
column. Then, the shift value{s) of a circulant is defined as the index (indices) of the nonzero 
element(s) in the first column. Note that a shift value takes the value from to 2; — 1 and 00 is 



used as a shift value of a zero matrix H, 



QC LDPC codes can be fully represented by binary polynomials as shown in IfTSl . This 
polynomial representation is based on the isomorphism between z x z binary circulants and the 
polynomial ring F2[a;]/(s^ + 1). The polynomial parity-check matrix H(x) of C is defined as 

H[x) = 

hj~ifi{x) hj_i^i{x) ■■■ /i,/_i,L-i(a;) 

where hj^i{x) = J2i=o ^ F2[a;]/(a;^ + 1) and hj^i^i is the element with the index i in the 

first column of Hj^i. We can see that the number of nonzero terms in hj^i(x), which is denoted 
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by wt{hj^i{x)), is equal to wt{Hj^i) and the degrees of all nonzero terms in hj^i{x) are equivalent 
to the shift values of Hj i. 

The protograph [111 of a QC LDPC code C is a bipartite graph whose incidence matrix is 
P = [pj^i], where pj^i = wt{Hj^i). There are two kinds of nodes in the protograph, where 
horizontal (check) nodes correspond to rows in P and vertical (variable) nodes correspond to 
columns in P. The Tanner graph of C is constructed by copying the protograph z times and 
cyclically permuting the same z edges. Such copy-and-permute operation is called lifting and 
the length of a subblock z is also called the lift size of C. If pj^i > 2, there are multiple edges 
between the horizontal node with index j and the vertical node with index I in the protograph. 
A shift value is assigned to each edge in the protograph so that an edge is lifted by using the 
cyclic permutation with the assigned shift value to generate C. Note that, in this paper, the term 
'protograph' refers to both the bipartite graph and its incidence matrix based on their equivalence. 

B. Inevitable Cycles 

Necessary and sufficient conditions on the existence of cycles in the Tanner graph of QC 
LDPC codes are derived in terms of shift values in These conditions are only applied to 
single-edge protographs but they can be naturally extended to cover the case of multiple-edge 
protographs as in Lemma [1} 

Let G = (V, E) denote a graph with a set of vertices V and a set of edges E. Let Vk (e^) 
represent a vertex (an edge) in V (E). A walk is an alternating sequence of vertices and edges, 
denoted by v^^e^^vi^ ■ ■ ■ Vi^_-^ei^_^Vi^, where the vertices Vi^ and vi,^^^ are the endpoints of the 
edge Ciy The length of a walk W, denoted by 1{W), is defined as the number of edges in W. 
A walk is closed if Vi^ = Vi^^ and a walk is non-reversing if e^^ 7^ ^i^+i f*^^ j = 0,1, . . . ,n — 2. 
A closed walk is said to be tailless if ej,^_^ 7^ e^y In this paper, connected graphs are only 
considered and it is noted that a cycle is defined as a closed walk whose traversed vertices and 
edges are all distinct and the length of the shortest cycle in a graph is called the girth of the 
graph. 

Cycles in the Tanner graph of a QC LDPC code are closely related to tailless non-reversing 
closed (TNC) walks in its protograph. The shift sum of a walk in a protograph, denoted by 
s{W), is defined as the alternating sum of shift values assigned to the edges in W, that is, s{W) = 
ES^"^(-1)^ (shift value of e,J. Lemmajl] shows necessary and sufficient conditions for a cycle 
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(a) An inevitable cycle of length 12. 
Fig. 1. Examples of inevitable cycles in QC LDPC codes. 



(b) An inevitable cycle of length 10. 



of a certain length in the Tanner graph of QC LDPC codes to be generated from the protographs, 
which can be applied to both single-edge protographs and multiple-edge protographs. Its proof 
is directly derived from the results in flU and [[TSl . 

Lemma 1: Let W denote the set of all TNC walks of length n in a protograph. Suppose that 
a QC LDPC code is lifted from the protograph with the lift size z. Then, the Tanner graph of 
this QC LDPC code has a cycle of length n if and only if there exists a walk W eW such that 
s{W) = mod z and W does not contain any shorter TNC walks with the zero shift sum. 

The girth of QC LDPC codes is determined by the structure of the protograph, the lift size, 
and all the shift values assigned to edges. However, we can derive an upper bound on the girth 
of QC LDPC codes lifted from protographs without considering the lift size and the shift values 
based on the concept of inevitable cycles [[3]|, flU, ifTOl . 

Definition 1: An inevitable cycle induced by a protograph is defined as the cycle which always 
appears in the QC LDPC code lifted from the protograph regardless of the lift size and the shift 
values. 

It is well known that a QC LDPC code whose protograph has the 2 x 3 (or 3 x 2) all-one 
matrix as its submatrix must have the inevitable cycles of length 12 [O, |j4l. In other words, the 
girth of this QC LDPC code is less than or equal to 12. Such an inevitable cycle of length 12 is 
depicted in Fig. |l(a) Also, in QC LDPC codes lifted from multiple-edge protographs, inevitable 



cycles can be induced. As an example. Fig. 1(b) shows an inevitable cycle of length 10, which 
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appears in QC LDPC codes lifted from protographs with double edges. We can see that for a 
certain subgraph structure, inevitable cycles are always generated no matter what shift values 
are assigned to edges. 

III. Subgraphs of Multiple-Edge Protographs Inducing Inevitable Cycles 

In order for QC LDPC codes to have large girth, their protographs should not contain the 
subgraphs which induce short inevitable cycles in the QC LDPC codes and thus it is necessary to 
find out all such subgraphs. From now on, the terms 'an inevitable-cycle-inducing (ICI) subgraph 
of length 2i' will refer to a subgraph inducing inevitable cycles of length 2i . In dH, ICI subgraphs 
of length up to 12 in single- and multiple-edge protographs were fully investigated and, in (TU\, 
all ICI subgraphs of lengths 12 to 20 in single-edge protographs were searched by a brute force 
method. After that, a graph-theoretical framework was provided in [fTSl . which can be used to 
search all single- and multiple-edge ICI subgraphs. In this section, we will search and provide 
all ICI subgraphs as an extension of flU, [fTOl . and [fTSll . 

Define V2i as a set of all irreducible ICI subgraphs of length 2i satisfying the following 
conditions: 

1) A subgraph P E V2i induces inevitable cycles of length 2i in the QC LDPC code. 

2) A subgraph P G V2i does not contain any proper subgraph which induces inevitable cycles 
of length less than or equal to 2i. 

3) The number of rows in a subgraph P G V2i is not larger than that of columns. 

4) From each isomorphic class in V2i, only one protograph must be chosen as a representative 
of that class. 

The conditions 1) and 2) guarantee that if a protograph does not have any subgraph P G V2i' 
for i' < i, the QC LDPC code appropriately lifted from this protograph has girth larger than 
or equal to 2i. A subgraph P G V2i takes an irreducible form because the condition 2) implies 
that if any edge is removed from P, it cannot induce inevitable cycles of length 2i. Conditions 
3) and 4) are required to choose a unique representative for each isomorphic class of subgraphs 
inducing inevitable cycles of length 2i. 

For identifying V2i, we need to investigate the relationship between inevitable cycles and 
TNC walks. A TNC walk W is called abelian-forcing [fTSl if for each edge in W, the number 
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(a) Theta graph. (b) Dumbbell graph. 

Fig. 2. Theta graph and dumbbell graph. 

of traversals of the edge in a direction is the same as that in the opposite direction. Clearly, the 
shift sum of abelian-forcing TNC walks is zero regardless of the shift values of their edges. An 
abelian-forcing TNC walk is said to be simple if it does not contain any shorter abelian-forcing 
TNC walks. It is obvious that inevitable cycles of QC LDPC codes are generated from simple 
abelian-forcing TNC (SAFTNC) walks in protographs. 

Lemma 2: Any abelian-forcing TNC walk contains at least two different cycles. 

Proof: Consider an abelian-forcing TNC walk W = Vi^Ci^Vi^ ■ ■ ■ Vi^_-^ei^_^Vi^. There exists a 
vertex Vj such that Vi^, = = vj for some k ^ I. Also, there exists a path Vmei^_^Vi^ ■ ■ ■ Vi^Ci^Vm 
in W such that all vertices from Vi^ to are distinct. Since W is non-reversing and tailless, 
that path forms a cycle and thus W contains at least one cycle. 

Assume that W contains only one cycle. Since W is abelian-forcing, there exists a path 
VfegVhei^_^Vi^ ■ ■ ■ Vi^eij^VhegVf in W such that Vi. ^ Vh for all j = a, a + 1, . . . ,b. This contradicts 
the assumption of W because W cannot move from a vertex to itself without reversing. Therefore, 
W contains at least two different cycles. ■ 

As in (15], two classes of graphs are defined as illustrated in Fig. |2] 

Definition 2 (/f75]/).- A (xi, X2, x^ytheta graph, denoted by T(xi, X2, x^), is a graph consisting 
of two vertices, each of degree three, that are connected to each other via three disjoint paths Xi, 
X2, X3 of the number of edges xi > 1, X2 > 1, and 0:3 > 1, respectively. A (zi, Z2; y)-dumbbell 
graph, denoted by D{zi, Z2;y), is a connected graph consisting of two edge-disjoint cycles Zi 
and Z2 of the number of edges zi > 1 and Z2 > I, respectively, that are connected by a path Y 
of the number of edges y > 0. 

Lemma 3: Connecting two different cycles always results in either a theta graph or a dumbbell 
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graph. 

Proof: Let Ci and C2 denote two different cycles. Then, Ci and C2 can be connected in 
only three ways: The number of common vertices in Ci and C2 is (i) 0, (ii) 1, or (iii) larger 
than or equal to 2. For the cases (i) and (ii), Ci and C2 form D{zi, Z2] y) with y > or y = 0, 
respectively. In the case (iii), T(xi, X2, ^3) is formed where Ci = Xi |J X2, C2 = X2 |J X3, and 
X2 + 1 is the number of the common vertices. ■ 

Lemma 4: The lengths of SAFTNC walks in T(xi, X2, X3) and D{zi, Z2] y) are 2(xi +X2 H-Xa) 
and 2(2:1 + Z2) + 4y, respectively. 



Proof: Consider T(xi, a;2, ^3) in Fig. 2(a) Let ui and U2 denote the left and the right 
vertices of degree three, respectively, and let Xi, X2, and X3 be the paths from ui to U2- Also, 
let Xi, X2, and X3 denote the reverse paths of Xi, X2, and X3, respectively. Then we can see 
that an SAFTNC walk X1X2X3X1X2X3 has the length 2{xi + X2 + X3) and any other SAFTNC 
walks possibly generated in T{xi, X2, X3) have the same length. 

Similarly, consider D(zi, Z2; y) in Fig. |2(b) Let vi and V2 denote the left and the right vertices 
of degree three, respectively, and let Zi and Z2 be the cycles rotating clockwise from vi and 
V2, respectively, and let Y be the path from vi to V2- Also, let Zi, Z2, and Y denote the reverse 
paths of Zi, Z2, and Y, respectively. Then we can see that an SAFTNC walk Z1YZ2YZ1YZ2Y 
has the length 2{zi + Z2) + 4?/ and any other SAFTNC walks possibly generated in -D(zi, Z2] y) 
have the same length. ■ 

Note that if any edge is removed from T(xi,X2,X3) or D{zi, Z2]y), those inherent SAFTNC 
walks disappear and thus T(xi, 0:2, ^3) and D(zi,Z2;y) are of irreducible form. Now, we will 
check whether it is sufficient to only consider theta graphs and dumbbell graphs for V2i- 

Lemma 5: Suppose that a graph G contains at least one of theta graphs or dumbbell graphs 
as its proper subgraphs. The shortest SAFTNC walk in G occurs only in a theta graph or a 
dumbbell graph. 

Proof: Let W denote the shortest SAFTNC walk and assume that W traverses all edges in 
G. From Lemmas |2] and [3| W should contain a theta graph or a dumbbell graph. Consider the 
following two cases: (i) G has some theta graphs, (ii) G does not have any theta graphs. 

In the case (i), we first note that l(W) is at least twice the number of edges in G due to the 
definition of abelian-forcing TNC walks. The SAFTNC walk only generated by a theta graph 
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in G is shorter than W because the SAFTNC walk has the length exactly twice the number of 
edges in the theta graph. This contradicts the assumption that W is the shortest one. In the case 
(ii), we note that a simple abelian-forcing TNC walk should traverse the edge not belonging to 
any cycles at least four times because if the walk traverses the edge twice, the walk will include 
two simple abelian-forcing TNC walks each of which occurs in the different side of the edge. 
Since cycles in G are connected with each other via only one path which does not belong to 
any cycles, the SAFTNC walk only generated by a dumbbell graph in G is shorter than W . 
This contradicts the assumption that W is the shortest one. Therefore, W occurs only in a theta 
graph or a dumbbell graph. ■ 

In the next theorem, V2i will be identified. 

Theorem 1: V2i is a collection of all T(xi, X2, a;3)'s with 2(xi + X2 + X3) = 2i and all 
D{zi, Z2] vYs with 2{zi + Z2) + Ay = 2i. 

Proof: From Lemmas [2] and [5| any subgraph P E V2i should be either a theta graph or a 
dumbbell graph. Therefore, the proof is completed by Lemma |4} ■ 

Now we can find all single- and multiple-edge ICI subgraphs from T{xi,X2, x^) and D{zi, Z2; y). 
A representative of an isomorphic class in V2i can be uniquely chosen by selecting parameters 
satisfying the following conditions: 

• Xi> X2> X^>\ 

• xi, X2, are all even or all odd 
. Zi>Z2>2,y>{) 

• zi and Z2 are even. 

Note that the second and the fourth conditions are derived because each subgraph P E V2i is a 
bipartite graph. 

According to Theorem [T| each integer solution of the equations 2{xi + X2 + ^3) = 2i and 
2(^1 + ^2) + % = 2i forms one ICI subgraph in V2i- Note that all ICI subgraphs of any length 
can be easily found and T{xi, 1, 1) and D{zi, 2; y) are ICI subgraphs having multiple edges. All 
ICI subgraphs of length up to 20 are listed as a form of theta or dumbbell graphs in Table |I] and 
all ICI subgraphs of length up to 14 are listed as a form of incidence matrices as follows: 

n = [3]; 
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TABLE I 

All ICI Subgraphs of Length up to 20 (T: Theta Graph, D: Dumbbell Graph, S: Single-Edge, M: 

Multiple-Edge) 



V2^ 


Pa 


Ps 


Pio 


Pl2 


Pii 


Pie, 


Pis 


P20 


T 


Xl 


1 




3 


2 




3 


5 


4 








3 


5 


7 


4 


6 










X2 


1 




1 


2 




3 


1 


2 








3 


3 


1 


4 


2 










X3 


1 




1 


2 




1 


1 


2 








3 


1 


1 


2 


2 










Type 


M 




M 


S 




S 


M 


S 








S 


S 


M 


S 


S 










D 


Zl 




2 




2 


4 






2 


4 


4 


6 








2 


4 


4 


6 


6 


8 


Z2 




2 




2 


2 






2 


2 


4 


2 








2 


2 


4 


2 


4 


2 


y 









1 









2 


1 














3 


2 


1 


1 








Type 




M 




M 


M 






M 


M 


S 


M 








M 


M 


S 


M 


S 


M 



2 2 



10 



12 



14 



2 1 

1 1 

2 1 

2 

2 1 

1 1 
1 1 



2 11 
1 1 



1 1 1 
1 1 1 



1 1 1 
1 1 
1 1 



where the single-edge ICI subgraphs in Table |I] were also listed in [[TOll and all ICI subgraphs 
of length up to 12 were also listed in [9|. Note that the transpose of each ICI subgraph also 
generates inevitable cycles of the same length and thus V2i will be used to denote both the listed 
matrices and their transposes. 
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IV. Construction of Regular Protographs Avoiding Inevitable Cycles of 

Length Less Than 12 

In this section, we will construct regular protographs which avoid inevitable cycles of length 
less than 12 in QC LDPC codes. Consider a regular J x L protograph of which the column- 
and row-weights are and dc, respectively, where J < L. If triple or more edges exist in the 
protograph, the girth of the lifted QC LDPC code is limited to 6 because of Vq = [3]. Therefore, 
only protographs with single and double edges will be considered in this paper. Let n2 denote 
the number of double edges in the protograph. 

Most of the considered protographs have at least two cycles and thus they always induce some 
inevitable cycles according to Lemmas |3] and |4j Note that even if a protograph is designed not 
to contain any 7^2?' with i' < i so that inevitable cycles of length less than 2i are avoided, this 
protograph may have some inevitable cycles of length larger than or equal to 2i. 

To construct protographs which do not induce inevitable cycles of length less than 10, a pair 
of 2's should not appear in any row or in any column of the protograph to avoid Vs- As in 
the next lemma, the number of double edges in a protograph should be upper bounded by the 
number of horizontal nodes to construct QC LDPC codes with girth larger than or equal to 10. 

Lemma 6: If a J x L protograph does not induce inevitable cycles of length less than 10, 
then n2 < J. 

Proof: If n2 > J, there always exists a row which has at least two 2's and thus the protograph 
contains Vs- This contracts the assumption. ■ 

In order for QC LDPC codes to have the girth larger than or equal to 12, their protographs 
should not contain Vq, Vs, and Viq. We will explain that an incidence matrix of a balanced 
ternary design (BTD) with p2 = 1 and A = 2 is also the incidence matrix of a regular protograph 
with n2 = J that does not induce inevitable cycles of length less than 12. 

Definition 3 ([19]): A balanced ternary design BTD(f , 6; pi, "^'i k, A) is an arrangement of 
V elements {1,2, ...,t>} into b multisets, or blocks, each of cardinality k, k < v, satisfying 
that (i) each element appears r = pi + 2p2 times altogether, with multiplicity one in exactly pi 
blocks, with multiplicity two in exactly p2 blocks and (ii) every pair of distinct elements appears 
A times, i.e., if rrij^h is the multiplicity of the element j in the h-th block, then for any elements 
i and j with i ^ j, we have X]/i=i ^i,hfnj,h = A. 
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TABLE n 

Regular Protographs With n2 = J Avoiding Inevitable Cycles of Length < 12 Constructed From BTDs 

FOR dc < 15 



J 


6 


12 


9 


20 


12 


30 


42 


48 


42 


15 


60 


L 


12 


24 


27 


40 


48 


60 


63 


64 


84 


75 


100 


dv 


3 


4 


3 


5 


3 


6 


8 


9 


7 


3 


9 


dc 


6 


8 


9 


10 


12 


12 


12 


12 


14 


15 


15 



Note that a v x b incidence matrix of a BTD(t>, b; pi, p2, r; k, A) is simply expressed as [nij^h] 
and the column- and row-weights are k and r, respectively. 

Theorem 2: An incidence matrix of a BTD{v,b; pi, p2,r;k, \) with p2 = I and A = 2 does 
not contain Vq, Vs, and Vio- 

Proof: Let Prtd be an incidence matrix of a BTD{v,b; pi, p2,r] k, X) with p2 = 1 and 
A = 2. Since every element of this BTD can have multiplicity up to two, Vq does not appear 
in Pbtd- The condition p2 = I implies that 2 appears once in each row of Pbtd and A = 2 
implies that each column of Pbtd can have at most one 2. Hence Pbtd does not contain Vs- 
Since a pair of distinct elements will appear at least three times in this BTD if Pbtd has Vio 
as its submatrix, Pbtd does not contain Vw- ■ 

All possible BTDs with r < 15 are given in [|20l . Table |Il] lists all parameters of regular 
protographs with dc < 15 avoiding inevitable cycles of length less than 12 constructed from 
BTDs. 



Example 1: An incidence matrix of BTD(6, 12; 4, 1, 6; 3, 2) is shown in Fig. 3(a) and we can 
see that any ICI subgraph V2i for i <h does not appear. 

As in Table |ll| the incidence matrices of BTDs with p2 = 1 and A = 2 do not provide 
a sufficiently large number of regular protographs. In fact, the condition that every pair of 
distinct elements appears exactly twice is not necessary and the condition that every pair of 
distinct elements appears at most twice is enough for constructing regular protographs avoiding 
inevitable cycles of length less than 12. Besides the regular protographs in Table |n| there are 
many regular protographs with n2 = J avoiding inevitable cycles of length less than 12. 
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(a) 6 X 12 (b) 6 X 8 

Fig. 3. Two regular protographs with dv = 3 and n2 = 6 avoiding inevitable cycles of length < 12. 



Example 2: Find the smallest regular protograph with = ?> and n2 = J avoiding inevitable 
cycles of length less than 12. We first derive a necessary condition for the existence of such a 
regular protograph by regarding the protograph as an incidence matrix of a block design. There 
are (2) distinct pairs of elements and, on the other hand, the number of all possible pairs of 
elements in the design is ^2 ■ 2 + (L — ^2) ■ (2) • Since every pair of elements appears at most 
twice, that is, Viq does not appear in the protograph, we have a necessary condition 

2- (2) >^2-2 + (L-n2)- Q- 

For J = 3, due to L > 4 and n2 = 3, the necessary condition Q is not satisfied. For J = 4, by 
counting the edges in the protograph, the equality Jd^, = d^L, that is, Adc = SL holds. Since the 
smallest integer root of this equality is {dc, L) = (6, 8), we have L > 8 and ([T]) is not satisfied. 
Similarly, for J = 5, L should be larger than or equal to 10 and ([T]) is not satisfied either. 

For J = 6, from 6dc = 3L, the possible smallest protograph has the size 6x8 and it satisfies 
([T]). By first constructing a 6 x 6 regular matrix where each column has one 2 and then properly 
adding two columns only consisting of O's and I's, a 6 x 8 regular protograph can be constructed 
as given in Fig. |3(b) This is the smallest regular protograph with d^ = S and n2 = J avoiding 
inevitable cycles of length less than 12. 



V. Construction of Regular Protographs Avoiding Inevitable Cycles of 

Length Less Than 14 

Now we will focus on the construction of regular multiple-edge protographs avoiding in- 
evitable cycles of length less than 14. A systematic construction method of single-edge regular 
protographs avoiding inevitable cycles of length less than 14 was provided in [ 10]. Since multiple- 
edge protographs are now being considered, two additional ICI subgraphs having double edges 
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Fig. 4. The structure of regular protographs avoiding inevitable cycle of length < 14. 

of Vi2 as well as and Pio must be avoided, which makes the problem more complicated. In 
this section, systematic construction methods for multiple-edge protographs are proposed based 
on various combinatorial designs. 

Consider a regular J x L protograph whose column- and row-weights are and dc, respec- 
tively. Let n2 denote the number of double edges in the protograph. Assume that dy > 3 because 
the regular QC LDPC codes with d^ = 2 are not used in general due to their poor performance. 

Using row and column permutations, every regular protograph not inducing inevitable cycles 
of length less than 14 can be represented as in Fig. |4} The n2 x n2 submatrix A has n2 2's as 
its diagonal elements and the other elements of A should be zero to avoid the first ICI subgraph 
of Vi2- F is a (J — 722) X submatrix consisting of columns of weight dy — 2. By appropriate 
column permutation of all but A and F in the protograph, all the columns whose lower parts 
have nonzero weight are relocated in the part of B and G, and the remaining columns make T 
with column-weight dy and all-zero matrix O. Let G and T be x Lq and Jt x Lt matrices, 
respectively. 

By Lemma |6} 722 cannot be larger than J. Moreover, if the regular protographs which do not 
induce inevitable cycles of length less than 14 are considered for d^ = 3, the following theorem 
provides an additional condition on n2. 

Theorem 3: Assume that a regular protograph with d^ = 3 and (ic > 4 does not induce 
inevitable cycles of length less than 14. Then n2 < J — 2. 

Proof: The inequality n2 < J holds by Lemma |6j The protograph with n2 = J should be 
of the form [yl|i?|r] from Fig. |4] and the submatrix A is no longer a diagonal matrix due to 
dv = 3. Therefore, A should contain the first ICI subgraph of V12 because each column of A also 
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contains exactly one 1 and hence n2 should be less than J. There should be one double edge in 
each row to avoid Vs and each column should have at most one double edge. The column with 
a double edge also has a single edge in other position, which generates the first pattern of Vu- 
Now suppose that n2 = J —1. The protograph has the form of Fig. |4] and F is the 1 x ( J — 1) 
all-1 matrix. Due to F, dc cannot be less than J— 1. If dc > J— I, G becomes the 1 x (dc — {J—l)) 
all-1 matrix and each column of B has a pair of I's, which generates Vio in the union of A, B, 
F, and G. If d^ = J — I, the protograph is made up of only A, T, F, and O, and the size of 

Tis (J— 1) X (J — 1)(J — 3)/3 because the column- and row- weights of T are 3 and J — 3, 

1 1 

respectively. Since T should not have as its submatrix to avoid the second and the third 

1 1 

ICI subgraphs of Vu in the union of A, B, and T, a pair of I's in the same column can appear 
at most once in T. To satisfy this condition, the number of all possible column- wise pairs of I's 
should be larger than or equal to the number of actual column-wise pairs of I's in T. Therefore, 
we have (-^2^) ^ (2) ' ~ ^)(^ " 3)/3, i.e., J < 4. Due to 4 = J - 1, this contradicts the 
assumption of dc> 4:. ■ 



Based on Theorem (Sj the case of d^ = S and n2 = J — 2 is considered in Subsection |V-A 



and the construction method of regular protographs for d^ = 3 and n2 = J — 2 is extended not 



only to the case of = 3 and n2 < J — 2 but also to the case of > 4 in Subsection V-B 



A. Regular Protographs With (it, = 3 and n2 = J — 2 

In this subsection, the construction of regular protographs with d^ = 3 and n2 = J — 2 
is elaborated. Necessary conditions on dc and J for the existence of regular protographs with 
d^ = 3 and 77,2 = J — 2, which avoid inevitable cycles of length less than 14, are derived as 
follows. 

Theorem 4: Assume that a regular protograph with d^ = 3, dc> 4, and n2 = J — 2 does not 
induce inevitable cycles of length less than 14. Then dc and J should satisfy either 



1) 


J = 


5 


mod 6, J > 11, and dc 


= (J + l)/2 or 


2) 


J = 


2 


mod 6, J > 14, and dc 


= (J - 2)/2 or 


3) 


J = 


3 


mod 6, J > 9, and dc = 


= (J - l)/2, (J+ l)/2 or 


4) 


J = 


1 


mod 6, J > 13, and dc 


= (J - l)/2 or 


5) 


J = 





mod 6, J > 12, and dc 


= (J - 2)/2, J/2 or 
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6) J = 10 and 4 = 6. 

Proof: By counting the edges in the protograph, we have df-J = d^L. Since dy = 3 and L 
is an integer, dc-J = mod 3. Also, the submatrix F in Fig.|4]is a 2 x (J — 2) matrix consisting 
of weight- 1 columns. Consider two cases: (i) F contains an all-1 row, (ii) F does not contain 
any all-1 rows. 

For the case (i), if dc > J — 2, the ICI subgraph Vio appears in the union of A, B, F, and 
G. If dc = J — 2, G is a 2 X (J — 2) matrix with an all-1 row at the different row position from 
the all-1 row of F. Then there exist some rows containing a pair of I's in i? because B has 
2(J — 2) I's and the column- weight of B is 2, which generates the second ICI subgraph of Vu 
in the union of A, B, and G. Therefore, the case (i) is impossible. 

For the case (ii), if a column of G has a pair of I's, the column including this pair in the 
protograph and another column in the union of A and F generate Vw- Therefore, each column 
of F and G cannot have a pair of I's. Since the number of columns in F is J — 2 and the total 
number of columns in F and G is 2dc, we have (J — 2)/2 < d^., where the equality holds when 
B and G do not appear in the protograph. If a row of B has a pair of I's, either Vio or the 
second ICI subgraph of Vu must occur in the union of A, F, G, and B. Therefore, each row of 
B can have at most one 1 so that the number of I's in i? cannot exceed the number of rows in 
B. Since the column-weight of is 2 and B has 2(2^- (J-2)) I's, we obtain 4 < 3(J-2)/4 

from 2(2dc — (J — 2)) < J — 2. Finally, it remains to determine the structure of T such that the 

1 1 

submatrix does not appear in the union of B and T to prevent the second ICI subgraph 

1 1 

of Vi2- As m the proof of Theorem [3} by counting the number of column-wise pairs of I's in 
B and T, we obtain the condition Q ■ (24 - (J - 2)) + Q ■ {dJ/S - 24) < (^2^) yielding 
4< (J-l)(J-2)/(2(J-4)). 

The above conditions on dc and J are summarized as follows: 

dcJ = mod 3; 

J -2 , (3,^ , (J-l)(J-2) 

< 4 < min<^ -(J-2)/ ^ 



2 - [A^ " 2(J-4) 

Since dc and J are integers, the above conditions reduce to simple linear relations with respect 
to J modulo 6 as given in the theorem statement. ■ 
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TABLE III 

All Possible Regular Protographs Avoiding Inevitable Cycles of Length < 14 When di, = 3 and n2 = J — 2 

FOR J < 26 



J 


9 


10 


11 


12 


13 


14 


15 


16 


17 


L 


12, 15 


20 


22 


20, 24 


26 


28 


35, 40 




51 


dc 


4, 5 


6 


6 


5, 6 


6 


6 


7, 8 




9 


J 


18 


19 


20 


21 


22 


23 


24 


25 


26 


L 


48, 54 


57 


60 


70, 77 




92 


88, 96 


100 


104 


dc 


8, 9 


9 


9 


10, 11 




12 


11, 12 


12 


12 



In Theorem |4} all possible regular protographs avoiding inevitable cycles of length less than 



14 are provided for = 3 and ^2 = J — 2, and Table III only lists those for J < 26 among 
them. 

Now we focus on the existence problem and the construction of the regular protographs with 
the parameters found in Theorem |4j Note that the proposed protographs we will construct may 
not be all instances with the parameters in Theorem |4] but we provide at least one instance per 
each set of parameters and also note that Jg = 2, Lq = 2dc — {J — 2), Jt = J — 2, and 
Lt = dc{J — 6)/3. For given J and d^, the matrices B, T, F, and G can be constructed step 
by step as follows: 

1. For constructing B and T at once, an incidence matrix of a combinatorial block design 
suitably chosen for each case in Theorem |4] is modified such that it has the size Jt x {Lq + 
Lt), each of Lq columns corresponding to B has a disjoint pair of I's, the other columns 
have the weight 3, all rows have the weight dc — 2, and any column- wise pair of I's appears 
at most once to avoid the second and the third ICI subgraphs of V12 in [y4|i?|T]. 

2. In G, I's are placed such that [Lg/2\ columns have I's in the first row and the other 
columns have I's in the second row. 

3. For constructing F, I's are placed such that the union of A, B, F, and G does not contain 
Vw 

Note that the placement of I's in the third step is guaranteed by the bound dc < 3(J — 2)/4 in 
the proof of Theorem |4| 
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i.e., the subgraph 



. For this, the following Steiner system can be used. 



Since the conditions in Theorem |4] are necessary ones for the existence of T, a protograph 
may not exist for some parameter values. Therefore, we will show that there exist protographs 
for all parameter values given in Theorem |4] by providing explicit construction methods of B 
and T using various combinatorial designs as follows. 

1) J = 5 mod 6 and J > 11: 

In this case, we have dc = ( J + l)/2, Lq = 3, and Lt = {J + 1){J — 6)/6. We need to 
construct [B\T] of size (J — 2) x (J^ — 5J + 12)/6 to avoid repeated column-wise pairs of I's, 
1 1 
1 1 

Definition 4 (^T9\l): A t-{v, k, A) design is a pair {V, B), where is a i;-set of points and B 
is a collection of fc-subsets (blocks) of V with the property that every t-subset of V is contained 
in exactly A blocks in B. A Steiner system S{t, k, v) is the t-{v, k, A) design with A = 1. 

Lemma 7 f/fJPl/).' There exists S{2,3,v) only when v = 1,3 mod 6. 

The number of blocks in S{2,3,v) is v{v — l)/6. Since three columns have the weight two 
and the other columns have the weight three in the ( J — 2) x (J^ — 5 J + 12)/6 matrix 
the (J — 2) X ( J — 2)( J — 3)/6 incidence matrix of 5(2, 3, J — 2) may be modified to be used 
as [B\T] by deleting one 1 from each of well-chosen three columns and adding one column of 
weight three. In order for such a modified matrix to be a valid we should check whether 

three column- wise pairs of I's in the weight-2 columns are disjoint, all rows have the weight 
(J — 3)/2, and any column- wise pair of I's appears at most once. 

Without loss of generality, let {vi,V2,Vi}, {v2,V3,Vj}, and {fi,f3,ffc}, i ^ j ^ k, he 
three blocks of S'(2, 3, J — 2) corresponding to three columns containing a cycle of length 
6. Three disjoint blocks {v2,v.i}, {v3,Vj}, and {vi,Vk} are obtained by removing vi, V2, and 
f3 from {vi,V2,Vi}, {f2,f3,fj}, and {vi,V3,Vk}, respectively. Inserting a block {vi,V2,Vs} to 
this modified S{2, 3, J — 2) still makes every pair appear at most once. An incidence matrix 
of 5'(2, 3, J — 2) has the row-weight (J — 3)/2 and the above modification clearly keeps the 
row-weight unchanged. Therefore, we propose a construction method of [B\T] in the case of 
J = 5 mod 6 and J > 11 as follows: 

1. Permute the columns of an incidence matrix of S'(2, 3, J — 2) so that the first three columns 
contain a cycle of length 6. 
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2. Delete a 1 on the cycle of length 6 from each of the first three columns so that the resulting 
three column-wise pairs of I's are disjoint. 

3. Insert one column of weight three where three I's are located in the rows passed through 
by the above cycle of length 6. 

Actually, it is easy to choose three columns which contain a cycle of length 6 because an 
incidence matrix of 5'(2,3, J — 2) has many cycles of length 6. The following lemma shows 
how many cycles of length 6 exist in an incidence matrix of 5(2,3, J — 2). 

Lemma 8: An incidence matrix of 5(2, 3, J — 2) has ( J — 2) ( J — 3) ( J — 5)/6 cycles of length 

6. 

Proof: Consider three points vi, f 2, G ^ of 5(2, 3, J— 2). Three pairs {f 1, f 2}, {f 2, f 3}, {fa, f 1} 
appear in 5(2, 3, J— 2) in either of two ways: (i) one block has all the three pairs, that is, consists 
of f 1, V2, f 3; or (ii) each pair is contained in a block which does not have the other two pairs, that 
is, there are three blocks {f 1, f2, "Wj}, {f 2, v^^^Vj}, {^3, vi, Vk], where f j, Vj, Vk E V and i ^ j ^k. 
Three pairs in the case (ii) form a cycle of length 6 in the incidence matrix of 5(2, 3, J — 2). 
Hence the number of cycles of length 6 in the incidence matrix can be enumerated by substracting 
the number of all blocks from the number of the ways of choosing three points in V . This yields 
(V) - - 2)(^ - 3)/6 = (J - 2)(J - 3)(J - 5)/6. ■ 

Example 3: Fig. [5] illustrates the construction of an 11 x 22 protograph with = 3 and 
n2 = 9. A cycle of length 6 is denoted by the circles in the incidence matrix of 5(2, 3, 9) 



which has been already column-wisely permuted in Fig. 5(a) To obtain [-B|T], three I's marked 
by dotted circles are deleted and the column with I's in the first, the second, and the fourth 
rows is inserted as the first column of T. Let Vi, i = 1, . . . ,9, denote the points of 5(2, 3, 9), 
which also denotes the i-th row of [-B|T]. We can see that [B\T] does not have three pairs of 
I's {vi,V3}, {v4,v-/}, {v2,Vg} in any column and three pairs {^2,^3}, {vi,Vj}, {v4,Vg} in B are 
disjoint. The resulting 11 x 22 protograph with dy = 3 and 722 = 9 is shown in Fig. 5(b)| and we 
can check that V2i with i < 6 does not appear in this protograph. 

2) J = 2 mod 6 and J > 14: 

In this case, we have d^ = (J — 2)/2, Lq = 0, and Lt = (J — 2)(J — 6)/6. Since B and G do 
not appear in the protograph, T should be designed to avoid repeated column-wise pairs, where 
T has constant row- weight (J — 6)/2 and column-weight 3. A configuration whose incidence 
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(b) An 11 X 22 regular protograph with — 3 and n2 — 9. 
Fig. 5. The construction of an 11 x 22 regular protograph with dv — 3 and n2 — 9. 



matrix has the column-weight 3 and the size (J — 2) x (J — 2)(J — 6)/6 can be used for T. 

Definition 5 ([19]): A configuration (vr, h^) is an incidence structure of v points and h blocks 
such that (i) each block contains k points, (ii) each point lies on r blocks, and (iii) two different 
points are contained in at most one block. \i v = h and hence r = k, the configuration is called 
symmetric and denoted by Vk- 

It is important to check the existence of the configuration with the required parameters. 
The following theorem shows that such configuration always exists and therefore T can be 
constructed. 

Theorem 5: There exists a configuration {vr,bk) with v = J — 2, b = ( J — 2)( J — 6)/6, = 3, 
and r = ( J - 6)/2 for all J = 2 mod 6 and J > 14. 

Proof: Necessary conditions for the existence of {vr,bk) configuration [21] are given as (i) 
V <b and k < r, (ii) vr = bk, and (iii) v > r{k — l) + 1. We can easily check that the parameters 
in the theorem statement satisfy these conditions. Finally, the existence of such configurations is 
guaranteed by Theorem 3.1 in lillL that is, there exists a configuration with A; = 3 if and only 
if the necessary conditions hold. ■ 



October 16, 2012 



DRAFT 



22 



Now, a construction method of T is proposed based on the results in [llTTl . which uses 
configurations with parallel classes and resolvable configurations. 

Definition 6 (fd9^): A parallel class in a design is a set of blocks that partition the point set. 
A resolvable design is a design whose blocks can be partitioned into parallel classes. 

For J = 2 mod 6, S'(2, 3, J — 1) exists by Lemma [vj For J > 20, an incidence matrix 
of a resolvable configuration {vr,bk) with v = J — 2, b = (J — 2)(J — 4)/6, k = 3, and 
r = (J — 4)/2 can be constructed by removing a row and its incident columns in an incidence 
matrix of S{2, 3, J — 1) [21]. For J = 14, there is no resolvable configuration (I25, 2O3) but we 



can find a configuration (I25, 2O3) in the same manner as illustrated in Fig. 6(b) which contains 
some parallel classes from S{2, 3, 13) [|2T|. Since a parallel class of a configuration {vr, bk) with 
v = J-2, b={J- 2)(J - 4)/6, A; = 3, and r = (J - 4)/2 consists of (J - 2)/3 blocks 
and has all points exactly once, we obtain T by removing one parallel class from the incidence 
matrices of these configurations. The construction procedure of T for J = 2 mod 6 and J > 14 
is summarized as: 

1. Construct 5(2,3, J - 1). 

2. Make an incidence matrix of a resolvable configuration {vr,bk) with v = J — 2, b = 
(J — 2)(J — 4)/6, k = 3, and r = (J — 4)/2 by removing a row and its incident columns 
in an incidence matrix of S{2, 3, J — 1). 

3. Remove one parallel class which consists of (J — 2)/3 columns to obtain T. 



Example 4: An incidence matrix of 5(2, 3, 13) is shown in Fig. 6(a) An incidence matrix of 



a configuration (I25, 2O3) in Fig. 6(b) is constructed by removing the eighth row and its incident 



columns in the incidence matrix of 5(2,3,13) in Fig. 6(a) We see that the fourth, the sixth, 
the thirteenth, and the sixteenth columns form a parallel class. By removing these columns, an 
incidence matrix of a configuration (I24, 163) is constructed, which is used as T. The resulting 



14 X 28 protograph with d.^ = 3 and = 12 is shown in Fig. 6(c) 

3) J = 3 mod 6 and J > 9: 

3.1) dc = {J - l)/2 except for J = 9; 

In this case, we have Lq = 1 and Lt = (J — 1){J — 6)/6, and thus B should have only one 
pair of I's. Since 5(2, 3, J - 2) exists by Lemma |7| [B\T] may be constructed by removing 
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(a) An incidence matrix of S(2, 3, 13). 
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(b) An incidence matrix of a configuration (125,203). 
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(c) A 14 X 28 regular protograph with d„ = 3 and n2 — 12. 
Fig. 6. The construction of a 14 x 28 regular protograph with d„ = 3 and n2 = 12. 



J/3 — 1 columns from a (J — 2) x (J — 2)(J — 3)/6 incidence matrix of S{2, 3, J — 2) and 
then deleting a 1 in other column. To achieve the desired row-weight (J — 5)/2 of [i?|T], there 
should exist a submatrix consisting of J/3 columns satisfying that (i) two rows have weight 2 
and the others have weight 1 and (ii) one column has 1 at each of these two rows of weight 2. 

As seen in the case of J = 2 mod 6 and J > 14, for J = 3 mod 6 and J > 15, S{2, 3, J— 2) 
contains as a substructure a configuration (vr, bk) with v = J — 3, b = ( J — 3)( J — 5)/6, k = 3, 
and r = (J — 5)/2 which has at least one parallel class consisting of J/3 — 1 blocks. This implies 
that there are J/3 — 1 blocks in 5(2,3, J — 2), which partition all but one point. Also, there 
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(a) An incidence matrix of S{2, 3, 13). 
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(b) A 15 X 35 regular protograph with d^j — 3 and n2 — 13. 
Fig. 7. The construction of a 15 x 35 regular protograph with d„ = 3 and n2 — 13. 



always exists another block containing that point in S{2,3,J — 2). These J/3 blocks satisfy 
the above requirements (i) and (ii) for [i?|T]. The construction procedure of [B\T] for J = 3 
mod 6, J > 15, and dc= (J — l)/2 is summarized as: 

1. Construct 5(2,3, J -2). 

2. Select one row in an incidence matrix of S* (2, 3, J — 2) such that if the row and its incident 
columns are removed from the incidence matrix, the remaining part forms an incidence 
matrix of a configuration {vr,bk) with v = J ~ 3, b = (J — 3)(J — 5)/6, A; = 3, and 
r = (J — 5)/2 including at least one parallel class. 

3. Find J/3 — 1 columns which form one parallel class in the above configuration. 

4. In the incidence matrix of S'(2, 3, J — 2), delete a 1 in the selected row in Step 2 and remove 
the J/3 — 1 columns obtained in Step 3. 

5. Move the column which had the deleted 1 in Step 4 to the leftmost to obtain 

Note that the above construction method cannot be applied to the case of J = 9 and (ic = 4 
because the configuration (62, 43) does not have any parallel class. The case of J = 9 and (ic = 4 
will be covered in the last part of this subsection. 
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Fig. 8. A 9 X 15 regular protograph with d„ = 3 and 712 = 7. 



Example 5: Fig. |7] illustrates the construction of a 15 x 35 regular protograph with = 3 and 
n2 = 13. In an incidence matrix of S{2, 3, 13) in Fig. 7(a) the fifth, the seventh, the sixteenth, 
and the twentieth columns partition the set of row indices except for the index of the eighth 
row and the fourth column has a 1 in the eighth row. Thus, the 1 in the fourth column and the 
eighth row is deleted and the four boxed columns are removed from the incidence matrix. Then 
the resulting column of weight 2 is moved to the leftmost and a 15 x 35 protograph with dy = 3 



and r?,2 = 13 is shown in Fig. 7(b) 
3.2) d, = {J+l)/2- 

In this case, we have Lq = 3 and Lt = (J + 1)(J — 6)/6, and there are three column-wise 
pairs of I's in B. Since S'(2, 3, J — 2) exists for J = 3 mod 6 and J > 9 by Lemma |7} the 
construction method for J = 5 mod 6 and J > 11 can also be applied to this case in the same 
way. As an example, a 9 x 15 regular protograph with d^ = 3 and ?22 = 7 is shown in Fig. [8j 

4) J = 1 mod 6 and J > 13: 

In this case, we have dc = {J — l)/2, Lq = 1, and Lt = {J — 1){J — 6)/6. To construct 
[i?|T], start with 5(2,3,7) which always exists by Lemma [t] Similar to the case of J = 2 
mod 6 and J > 14, a configuration (vr, hk) with v = J~l, h={J— 1)(J — 3)/6, k = 3, and 
r = (J — 3)/2 can be constructed by removing a row and its incident columns in an incidence 
matrix of 5(2,3, J). Since any (J — 3)/2 blocks sharing a common point partition all points 
except the common point and another point in the configuration, by removing any row and its 
incident columns in an incidence matrix of the configuration, a (J — 2) x (J — 3)(J — 4)/6 matrix 
with J — 3 rows of weight (J — 5)/2 and a row of weight (J — 3)/2 is obtained. Removing a 1 
in the row of weight (J — 3)/2 results in a matrix which has the desired row-weight (J — 5)/2 
and exactly one column of weight 2. Clearly, this matrix can be used as [-B|T]. The construction 
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(a) An 11 X 15 modified matrix. 
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(b) A 13 X 26 regular protograph with d^j — 3 and n2 = 11. 



Fig. 9. The construction of a 13 x 26 regular protograph with ci„ = 3 and n2 — 11. 



procedure of [B\T] for J = 1 mod 6 and J > 13 is summarized as: 

1. Construct a configuration {vr,bk) with v = J — 1, b={J— 1)(J — 3)/6, A; = 3, and 
r = ( J — 3)/2 from S{2, 3, J) similar to the case of J = 2 mod 6 and J > 14. 

2. Obtain a matrix of size (J — 2) x (J — 3)(J — 4)/6 by removing a row and its incident 
columns. 

3. Delete a 1 in the row of weight (J — 3)/2 and move the column having the deleted 1 to 
the leftmost to obtain [B\T]. 



Example 6: An incidence matrix of a configuration (125,203) is given in Fig. 6(b) By 
removing the first row and its incident columns, we obtain an 11 x 15 matrix shown in Fig. 9(a) [ 
where the seventh row has the weight 5 and the others have the weight 4. Then the 1 in the 
seventh row is deleted from the second column and the second column is moved to the leftmost. 



The resulting 13 x 26 protograph with = 3 and 722 = 11 is shown in Fig. 9(b) 

5) J = mod 6 and J > 12: 
5.1) 4= (J-2)/2; 
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In this case, we have Lq = and Lt = (J — 2)(J — 6)/6. Similar to the case of J = 2 mod 6 
and J > 14, an incidence matrix of a configuration (vr, hk) with v = J — 2, b = ( J — 2)( J — 6)/6, 
k = 3, and r = (J — 6)/2 can be used as T. Such configuration can be constructed by using 
difference triangle set (DTS). 

Definition 7 ([19]): An (n, m)-dijference triangle set, or (n, m)-DTS, is a set W = {f/i, . . . , f/„}, 
where fox 1 < i < n, Ui = {ajo,aji, • • • iCiim} with a^j an integer satisfying = < an < 
■ ■ ■ < aim, and the differences an — aij over the integers for all i, j, I, 1 < i < n, < I ^ j < m, 
are all distinct and nonzero. 

Theorem 6 f/EH/j-' If there is an [n, 2)-DTS, a configuration [vr, hk) for t> > 6n + 3, 6 = nv, 
k = 3, and r = nk can be constructed from this DTS. 

For J = mod 6 and J > 12, a configuration (t>.,., bk) with v = J — 2, b = ( J — 2) ( J — 6)/6, 
k = 3, and r = (J — 6)/2 can be constructed from ((J — 6)/6, 2) -DTS by Theorem |6] According 
to [1211 . the construction procedure of T is provided as: 

1. Construct a ((J — 6)/6, 2)-DTS with Ui = {aiQ, an, 042}, i = 1, . . . , (J — 6)/6. 

2. For each Ui, i = 1, . . . , {J — 6) / 6, construct a column of length J — 2 denoted by Ci, which 
has 1 at the (ojo + l)-st, the {an + l)-st, and the (ai2 + l)-st rows and at other rows. 

3. For each i, construct a (J — 2) x (J — 2) matrix whose j-th column, j = 1, . . . , J — 2, is 
obtained by cyclically shifting Ci downward j — 1 times. 

4. Concatenate (J — 6)/6 matrices in Step 3 to obtain T. 

Note that we can easily construct a ((J - 6)/6,2)-DTS from the lists of DTS in [|22l for 



J = mod 6 and J > 12. Fig. 10 shows a 12 x 20 regular protograph with = 3 and ^2 = 10 
constructed from (1,2)-DTS. 

5.2) 4 = J/2; 

In this case, we have Lq = 2 and Lt = J {J — 6)/6. Pairwise balanced designs (PBDs) can 
be used to construct [i?|T]. 

Definition 8 (lU^): Let K be a subset of positive integers and let A be a positive integer. A 
pairwise balanced design of order v with block sizes from K, denoted by FBD(v,K; A), is a 
pair (V,i3), where V is a point set of cardinality v and B is a family of blocks of V which 
satisfy that (i) if B E B, then \B\ E K and (ii) every pair of distinct elements of V occurs in 
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Fig. 10. A 12 X 20 regular protograph with dv — 3 and n2 — 10. 

exactly A blocks of B. 

Let FBD{v,K) denote a FBD{v,K;\) with A = 1 and use FBD{v,KU {k*}) to denote a 
PBD containing only one block of size k in the PBD, where k ^ K is a positive integer. For 
J = mod 6, it was shown in [23] that PBD( J — 1, {3, 5*}) always exists. Note that five rows 
sharing 1 with the column of weight 5 have the weight (J — 4)/2 and the other rows have the 
weight (J - 2)/2 in a (J - 1) x (J^ - 3 J - 12)/6 incidence matrix of PBD(J - 1, {3, 5*}). 

Theorem 7: Removing a row of weight ( J— 4)/2 and its incident columns except the weight-5 
column from an incidence matrix of PBD(J — 1, {3, 5*}) makes a ( J — 2) x (J^ — 6 J + 6)/6 
matrix of constant row- weight (J — 4)/2. 

Proof: Without loss of generality, assume that the first column has I's at the first five rows 
in an incidence matrix of PBD(J — 1, {3,5*}). Consider the (J — 1) x (J — 4)/2 submatrix 
which consists of the columns incident to the first row. Except the first row, each row of this 
submatrix has only one 1 because every column- wise pair of I's should appear exactly once in 
an incidence matrix of this PBD. Thus, the ( J — 2) x ( J — 6)/2 submatrix obtained by removing 
the first row and the first column from the (J — 1) x (J — 4)/2 submatrix does not have 1 in the 
first four rows and each of the other rows has only one 1 . After removing the first row and the 
(J — 2) X (J — 6)/2 submatrix from the incidence matrix of PBD( J — 1, {3, 5*}), the remainder 
forms the ( J - 2) x ( - 6 J + 6) /6 matrix of row-weight (J - 4) /2. ■ 

The matrix constructed in Theorem [v] cannot be directly used as [B\T] due to the improper 
number of columns and the weight-4 column, but it can be easily modified to meet the require- 
ments for [B\T] by splitting the weight-4 column into two weight-2 columns. The construction 
procedure of [B\T] for J = mod 6, J > 12, and dc = J/2 is summarized as: 
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(a) An incidence matrix of PBD(11, {3, 5*}). 
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(b) A 12 X 24 regular protograph with d^j — 3 and n2 — 10. 



Fig. 11. The construction of a 12 x 24 regular protograph with d„ = 3 and n2 — 10. 



1. Construct PBD(J - 1, {3, 5*}). 

2. Remove a row of weight (J — 4)/2 and its incident columns except the weight- 5 column 
from an incidence matrix of PBD(J — 1, {3, 5*}). 

3. Split the weight-4 column into two weight-2 columns and move them to the leftmost to 
obtain [B\T]. 

Example 7: The construction process for a 12 x 24 regular protograph with = 3 and n2 = 10 



is illustrated in Fig. 11 An incidence matrix of PBD(11, {3, 5*}) is shown in Fig. 11(a) We 
can see that the submatrix consisting of the columns incident to the first row has exactly one 1 
in each row except the first row. By removing the first row and the second, the third, and the 
fourth columns and splitting the weight-4 column into two weight-2 column, [B\T] is obtained. 



The resulting 12 x 24 regular protograph with dy = 3 and n2 = 10 is shown in Fig. 11(b) 

J = 9, 4 = 4 and J = 10, 4 = 6: 

There only remain two cases to provide the construction methods of all regular protographs 
in Theorem |4j When J = 9 and dc = 4, we have Lq = 1 and Lt = 4, and [E|T] is a 7 x 5 
matrix with row-weight 2. Although the construction method of [B\T] for J = 3 mod 6 and 
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(a) An incidence matrix of S{2, 3, 7). 
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(b) A 9 X 12 regular protograph with dv — 3 and n2 — 7. 



Fig. 12. The construction of a 9 x 12 regular protograph with d„ = 3 and 712 = 7. 



dc = {J — l)/2 cannot be directly used, we can construct [B\T] from an incidence matrix of 



S{2, 3, 7) in Fig. 12(a) Since any two columns of an incidence matrix of S{2, 3, 7) have a 
common 1, removing the first two columns from an incidence matrix results in a 7 x 5 matrix 
where one row has the weight 1, four rows have the weight 2, and the remaining two rows have 



the weight 3 as shown in Fig. 12(a) To obtain [-B|T], first delete 1 from each of two rows of 
weight 3 in the 7x5 matrix such that two deleted I's do not belong to the same column and 
the columns containing two deleted I's do not have 1 in the row of weight 1. These two deleted 



I's are marked by circle in Fig. 12(a) Then, by replacing a at the row of weight 1 and one of 
the columns containing the deleted I's with a 1, [B\T] is constructed and the resulting 9 x 12 



regular protograph is shown in Fig. 12(b) 



When J = 10 and dc = 6, we have Lg = 4 and Lt = 8, and B has disjoint four column-wise 
pairs of I's. An incidence matrix of a symmetric configuration 83 can be used as T, which does 
not have disjoint four column- wise pairs of I's. A 10 x 20 regular protographs with d^ = 3 and 
n2 = 8 is shown in Fig. [13] 



B. Regular Protographs With d^ = 3 and n2 < J — 2, and With d^ > 4 

Regular protographs which do not induce inevitable cycles of length less than 14 for the case 
of d^ = 3 and ^2 < J — 2 and the case of d^ > A also have the same structure in Fig. |4] 
and the construction method in the previous subsection can be similarly applied to these cases. 
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Fig. 13. A 10 X 20 regular protograph with dy — 3 and n2 — 8. 

However, we do not elaborate on deriving necessary conditions like Theorem |4] and providing 
specific construction methods for all cases because they should be done case by case and are very 
lengthy. Instead, for given J, L, dy, dc, and n2, we provide a general framework for checking 
the constructibility and constructing each submatrix. 

First, some basic conditions on the parameters J, L, d^, dc, and n2 are provided to determine 
whether a regular protograph with the given parameters can be potentially constructed. In F, 
the number of all possible column- wise pairs of I's should be larger than or equal to the 
number of actual column- wise pairs of I's, that is, ("'^"^) > n2(f"^'^). Also, the last J — n2 
rows must have (J — n2)dc I's and the matrix F must have n2{dy — 2) I's, and thus we have 
(J - ^2)4 > n2{dy - 2). 

Second, consider constructing [B\T] of size r?,2 x (L — ^2). The matrix [B\T] has the constant 
row- weight dc — 2 and does not have any repeated column- wise pairs of I's to avoid the second 
and the third ICI subgraphs of V12 in [A|i?|T], and T has the constant column- weight dy. The 
matrix [B\T] can be constructed from an incidence matrix of block designs such as 5'(2, k,v), 
configurations (vr,bk), FBD(v, K), group divisible designs (GDD) [19| and so on because they 
do not have any repeated column- wise pairs of I's. If an incidence matrix of an S{2, k,v) or a 
configuration {vr, b^) has the desired size of [i?|T], it can be directly used as [B\T] where there 
is actually no B. Otherwise, an incidence matrix of some block designs can be used as [B\T] by 
doing a simple modification. Note that an incidence matrix of PBDs or GDDs may not have a 
constant row-weight while an incidence matrix of 5(2, /c,f) or configurations (fr,&fc) is always 
regular. 

To obtain [B\T] from an incidence matrix of different size, we may use the following modi- 
fication schemes: 
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1) Remove some rows. 

2) Remove a row and some columns incident to the row. 

3) Remove some parallel classes. 

4) Delete some I's and insert some columns not to have any repeated column- wise pairs of 
I's. 

5) Insert some parallel classes not to have any repeated column-wise pairs of I's. 

By properly applying these modification schemes, an incidence matrix is changed into [B\T] 
having no repeated column-wise pairs of I's and proper size. Moreover, Schemes 1), 3), and 5) 
change all row-weights by the same amount and Schemes 2) and 4) flexibly control row-weights 
according to how to select columns and I's. Therefore, we can freely use the above modification 
schemes until the desired size and the constant row- weight — 2 of [B\T] are achieved. 

For given J, L, d^, dc, and n2, it may be possible for [B\T] to take various forms, which 
implies that each B may have a different number of columns and a different distribution of 
I's. Therefore, some bounds on Lq, or the number of columns in B, need to be derived to 
construct [B\T]. Since the number of I's in G is {J — n2)dc — n2{dy — 2) and each column 
in G can have weight from 1 to d^, we have Lq < — ^2)'^c — ^2('^d — 2) < dyLc which 
yields {{J — n2)dc — n2{dy — 2)}/dy < Lq < {J — n2)dc — n2{dy — 2). Therefore, when [B\T] is 
constructed by selecting and modifying an incidence matrix of a block design, the above bound 
on Lq must be considered. 

Lastly, consider constructing F of size ( J — ^2) x ^2 and G of size (7 — ^2) x Lq. Column- 
weights of G are already determined if B is designed and I's in G should be located to avoid 
the second and the third ICI subgraphs of Vu in the union of A, B, and G. Then, for a given G, 
{J — n2)dc I's in F should be located such that the union of A, F, and G does not contain the 
second and the third ICI subgraphs of V12, and the union of A, B, F, and G does not contain 
7^10 while enforcing row- weights of [F\G] to be dc and column- weights of F to be — 2. 

For given parameters J, L, d^, dc, and 712, a general procedure for constructing regular 
protographs which avoid inevitable cycles of length less than 14 is summarized as: 

1. Check if the parameters satisfy the conditions ("'"2"^) > ?^2(''Y^) ^^'^ {J — n2)dc > n2{dy — 
2). If the conditions are not satisfied, stop the procedure. 

2. Obtain Lq satisfying {{J — n2)dc — n2{dv — 2)}/ci„ < Lq < {J — n2)dc — n2{dv — 2). 

3. Construct [B\T] using Lq obtained in Step 2 from an incidence matrix of a proper block 
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Fig. 14. A 15 X 30 regular protograph with d„ = 3 and n2 — 12. 



design. 

4. Construct [F\G] satisfying the weight constraints such that the union of A, B, F, and G 
does not have Vio and V12 as its subgraph. 

Example 8: Consider the construction of a 15 x 30 regular protograph with dy = 3, dc = 6, 
and n2 = 12. The given parameters satisfy the conditions {'^~^^) > ?t^2(°'"2^^) and (J — n2)dc > 
^2(^1, — 2). The 12 X 18 matrix [B\T] has the row-weight 4 and 2 < Lq < 6. An incidence matrix 
of a symmetric configuration I23 is chosen for the construction of which is constructed 

by removing two parallel classes from a 12 x 20 incidence matrix of the configuration (I25, 2O3) 
in Fig. |6(b) By inserting a parallel class consisting of six weight-2 columns to an incidence 
matrix of the symmetric configuration I23, [B\T] with Lq = 6 is constructed, where any repeated 
column- wise pairs of I's do not appear. Since the column-weight of B is 2, the 3x6 matrix 
G should have the column- weight 1. Let each row of G have two I's. Then F should have the 
column- weight 1 and the row-weight 4, and the I's in F can be properly distributed so that Vio 
and V12 do not appear in the union of A, B, F, and G. The resulting 15 x 30 regular protograph 
with dv = 3, dc = 6, and n2 = 12 is shown in Fig. \\A\ 



Example 9: Consider the construction of a 28 x 49 regular protograph with dy = A, dc = 7, 
and 722 = 21. We can check that those parameters satisfy two necessary conditions for the 
construction and 2 < Lq < 7. For constructing a 21 x 28 an incidence matrix of a 

symmetric configuration 2I4 is considered. We can find a parallel class consisting of seven 
weight-3 columns such that if those seven columns are inserted to the incidence matrix, any 
repeated column- wise pairs of I's still do not appear. Thus, we obtain \B\T] with Lq = 7. The 
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Fig. 15. A 28 X 49 regular protograph with d„ = 4 and n2 — 21. 



7x7 matrix G should have the column-weight 1 and its row-weight can be set to 1 . The matrix 
F has the size 7 x 21, the column-weight 2, and the row-weight 6. Due to (2) = 21, F can be 
constructed not to have any repeated column-wise pairs of I's. Also, we can make F to avoid 
Vio in the union of A, B, F, and G. The resulting 28 x 49 regular protograph with = 4, 



dc = 7, and n2 = 21 is shown in Fig. 15 



VL Construction of QC LDPC Codes and Their Minimum Hamming Distances 
A. Construction of QC LDPC Codes From the Proposed Protographs 

To verify the effectiveness of the proposed protographs, QC LDPC codes will be constructed by 
determining the lift size and assigning an appropriate shift value to each edge of the protographs. 
Given a protograph, it is not easy to find all shift values even for a moderate lift size such that 
the girth of the QC LDPC code is the same as the length of the shortest inevitable cycle. Huang 
et al. [7] proposed a search algorithm for small lift size and a shift value assignment scheme 
to achieve the target girth based on greedy search. This algorithm is originally designed for 
single-edge protographs. However, by a slight modification, this algorithm can be extended to 
the case of multiple-edge protographs. 

Consider a J x L protograph P with the column-weight d^, the row-weight d^ and the lift size 
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z. Each column of P has shift values and let s/ j, i = 0, . . . , (i„ — 1, denote the i-th shift value 
of the /-th column in P. Our goal is to determine all shift values j and search the minimum 
z when a protograph and a target girth g of QC LDPC codes are given. Let W„ denote the set 
of all TNC walks of length n in P. Then, by Lemma [T| the condition for achieving the target 
girth g of QC LDPC codes is that for any W G Wn, n = 4, 6, . . . , - 2, the shift sum siyV) 
satisfies siyV) ^ mod z. However, it requires too much computational complexity to find si^i 
and the minimum z satisfying the above condition by considering all search space of ^ and z. 

In order to reduce the search space of shift values, let si^i = riimi as in 171, where is 
the {i + l)-st element of the set {0, 1, 3, 7, 12, 20, . . .} which is constructed from {ro = 0} by 
adding rj, i = 1, 2, . . ., in order such that = rj_i + min^ fc<j [N \ {\rj — rk\]]. Thus we only 
need to find L values of mi instead of d^L values of j. Moreover, for further reduction of 
computational complexity, mi is determined in a greedy manner, that is, shift values of the /-th 
column in P are determined by considering only the first / columns in P. For this, let Wn'' 
denote the set of all TNC walks of length n in the matrix consisting of the first / columns of P. 

For a given target girth g, if si^i is already determined such that s{W) ^ for any W G Wn, 
n = 4, 6, . . . , — 2, the minimum z, denoted by ^min, can be sub-optimally determined as 
Zmin = max I W G Wn, = 4, 6, . . . , (7 — 2} + 1. Note that for any z > Zmin, the target 
girth is achieved. 

An algorithm to construct QC LDPC codes of moderate length by determining all shift values 
and searching the minimum lift size, called Algorithm 1, is provided as follows. If the target 
girth g is set to the length of the shortest inevitable cycle, we can generate QC LDPC codes of 
moderate length with the maximum achievable girth from the proposed protographs. Note that 
the computational complexity of Algorithm 1 is the same for both single-edge protographs and 
multiple-edge protographs under the same parameter values. 

Four QC LDPC codes are generated by using Algorithm 1 . From the 9x15 protograph in Fig. [8} 
a (15000, 6000) QC LDPC code with girth 14, denoted by Proposed Code 1, is constructed, which 
has z = 1000 and {mj = {-105, 36, 45, 75, -69, -303, -393, 127, -31, -199, 200, 184, 86, 200, 



199}. From the 9x12 protograph in Fig. [T2(b)t a (3600, 900) QC LDPC code with girth 14, de- 
noted by Proposed Code 2, is constructed, which has z = 300 and {mi} = { — 12, 18, —39, 75, —57, 
120, 15, 17, 0, -6, -8, -8}. From the 6 x 12 protograph in Fig. |3(a)| a (7200, 3600) QC LDPC 
code with girth 12, denoted by Proposed Code 3, is constructed, which has z = 600 and 
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Algorithm 1: Greedy Search for the Minimum Lift Size and Shift Values 


INPUT: Target girth g, J x L protograph, search bound Fmax 




OUTPUT: (0 < / < L - 1) and z^r^in 




INITIALIZATION: tq = 0, = ri_i + mmj^k<i [N \ {\rj - rk\}] for 1 


< i < dfj ~ 1 


MAIN ROUTINE 




for / = to L — 1 begin 




for mi = -Finax to Tinax begin 




Let si^i = riUii for < i < — 1. 




If s{W) 7^ for any W E n = 4,6, . . . , g - 2, 




^SnK) =max||s(iy)| 1 W e n = 4,6, . . . , g - 


2) + 1. 


Otherwise, -2min("^0 = 




end 




Select the minimum -2min("^0 ^^'^ ^^^^ the minimum ^min("^0 


to z^^j^ and also save 


the argument to mi. 




If there are multiple minimums, randomly pick any one. 




end 




(L-l) 
^min ^min 





{mi} = {-93,7,47,-52,-29,-192,30,29,30,3,19,42}. From the 6 x 8 protograph in Fig. 
3(b)[ a (800,200) QC LDPC code with girth 12, denoted by Proposed Code 4, is constructed, 
which has z = 100 and {mi} = {-3, 85, -18, -6, -7, 2, -5, -15}. 

B. Upper Bounds on the Minimum Hamming Distance of the Proposed QC LDPC Codes 

Smarandache and Vontobel ifTSl derived two upper bounds on the minimum Hamming distance 
of QC LDPC codes. While one bound needs whole code specifications, e.g., the structure of the 
protograph, the lift size, and the shift values, the other bound only requires knowledge of the 
protograph. 

These two upper bounds are shown in Theorems [8] and [9| and they are directly derived by 
finding some low-weight codewords as in Lemma |9j Let Qs denote the submatrix of Q that 
contains only the columns of Q whose index appears in the set S. 
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Definition 9 f/lIU/j-' The permanent of an m x m matrix Q = [g^ j] over some commutative 
ring is defined to be 

perm(Q) := ^ JJ gi,<,(i) 

o- ie{0,...,m-l} 

where the summation is over all ml permutations a on the set {0, . . . , m — 1}. 

Lemma 9 (^W[l): Let C be a binary QC LDPC code defined by a J x L polynomial matrix 
H{x) with the lift size z. Let S be an arbitrary size-( J + 1) subset of {0, 1, . . . , L — 1} and let 

c{x) = [co(x), Ci{x), . . . , cl-i{x)], where Cj(x) is a polynomial over ¥2{x)/{x^ + 1) defined by 

perm , ii i e S 

Ci{x) = < 

I 0, otherwise. 

Then c(x) is a codeword of C. 

Theorem 8 (lU^}: Let C be a binary QC LDPC code defined by a J x L polynomial matrix 
H{x) with the lift size z. Then the minimum Hamming distance of C is upper bounded as 

dminiC) < _^^^min_^^ 2^^* (perm {Hs\{i}ix))) (2) 

nsi=j+i ies 

* 

where the operator min gives back the minimum value of all nonzero entries in a list of values. 

Theorem 9 (EIES): Let C be a binary QC LDPC code lifted from a J x L protograph P. Then 
the minimum Hamming distance of C is upper bounded as 

dmin{C) < min > perm (P5\{i}) . (3) 

SC{0,...,L-1} .^-^ ^ \IJ/ 

Theorems [8] and |9] imply that for given J, L, dy, and dc, these two upper bounds on the 
minimum Hamming distance of QC LDPC codes possibly increase as the number of multiple 
edges in the protograph increases, which is supported by examples for some regular protographs 
in ^TS^ . Note that the bound in ([3]) is not tighter than the bound in (|2]), but the former approaches 
the latter for a large z and proper shift values. 

Consider the (15000,6000) Proposed Code 1. The upper bounds in ([2]) and ([3]) for this code 
are 246 and 256, respectively. For comparison, a QC LDPC code with the same parameter values 
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001 100001001 



(a) 9 X 15 (b) 9 X 12 

Fig. 16. Two single-edge regular protographs with di, = 3 avoiding inevitable cycles of length < 14. 



is generated from the 9 x 15 single-edge regular protograph in Fig. 16(a) by using Algorithm 1. 
This single-edge protograph is constructed by attaching the last three columns to an incidence 
matrix of S'(2,3,9) to avoid inevitable cycles of length less than 14. The upper bounds in (j2]) 
and (|3]) for this code are 218 and 230, respectively. 

Consider the (3600,900) Proposed Code 2. The upper bounds in ([2]) and ([3]) for this code are 
362 and 416, respectively. For comparison, a QC LDPC code with the same parameter values is 
generated from the 9 x 12 single-edge regular protograph in Fig. 16(b)| by using Algorithm 1. By 
using the construction method in [fTOl . this single-edge protograph is constructed by concatenating 
an incidence matrix of a (92, 63) configuration and cyclically row-shifted matrix of it. The upper 
bounds in ([2]) and ([3]) for this code are 314 and 384, respectively. 

Consider the (7200,3600) Proposed Code 3. The upper bounds in ([2]) and ([3]) for this code 
are all 68. For comparison, a QC LDPC code with the same parameter values is generated from 
the 6 X 12 single-edge regular protograph in Fig. 17(a)| by using Algorithm 1. This single-edge 
protograph is the best one of randomly constructed protographs in the sense of upper bounds 
on the minimum Hamming distance. The upper bounds in ([2]) and ([3]) for this code are all 56. 

Finally, consider the (800,200) Proposed Code 4. The upper bounds in ([2]) and (|3]) for this code 
are 130 and 174, respectively. For comparison, a QC LDPC code with the same parameter values 



is generated from the 6x8 single-edge regular protograph in Fig. 17(b) by using Algorithm 1. 
This single-edge protograph is the best one of randomly constructed protographs in the sense 
of upper bounds on the minimum Hamming distance. The upper bounds in ([2]) and ([3]) for this 
code are 98 and 110, respectively. 

The above results clearly show that two upper bounds (|2]) and ([3]) on the minimum Hamming 
distance of QC LDPC codes are affected in a positive way by using double edges in the 
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101010101010 
010101100110 
100101010101 
010110010110 
011010011001 
101001101001 



1 1 1 1 
1 1 1 1 
10 110 1 
1 1 1 1 
1110 10 
10 10 10 1 



(a) 6 X 12 (b) 6 X 8 

Fig. 17. Two single-edge regular protographs with di, = 3 avoiding inevitable cycles of length < 12. 

protographs. In general, a multiple-edge protograph is more difficult to design than a single- 
edge protograph under the condition that they induce the shortest inevitable cycles of the same 
length. However, if multiple-edge protographs are once constructed, QC LDPC codes lifted from 
them can potentially give a larger upper bound on the minimum Hamming distance than those 
lifted from single-edge protographs. 



C. Comparison of Error Correcting Performance 

Performance of four proposed QC LDPC codes, that is. Proposed Code 1 to 4 is compared 
with those of the progressive edge-growth LDPC codes, called PEG 1 to 4 [|24| and the QC 
LDPC codes, called PEG QC 1 to 4 [|25l with the same code length, code rate, and column- 
weight. PEG LDPC codes and PEG QC LDPC codes are well known to have good error 
correcting performance comparable to those of random LDPC codes. Note that the girths of 
such (15000,6000), (3600,900), (7200,3600), and (800,200) PEG LDPC codes and PEG QC 
LDPC codes are 12, 12, 12, and 10, respectively, and these codes are obtained by the PEG 
algorithm to have as large girth as possible. 

Simulation is carried out through the binary input additive white Gaussian noise (BIAWGN) 
channel. The belief propagation (BP) decoding algoirthm is used and the number of maximum 
iterations is set to 100. The frame error rate (FER) performances of all the above LDPC codes 



are compared in Fig. 18 and we can see that the proposed QC LDPC codes show as good error 
correcting performance as the PEG LDPC codes and the PEG QC LDPC codes. Note that the bit 
error rate (BER) curves behave qualitatively the same as the FER curves and they are omitted 
in this paper. 



October 16, 2012 



DRAFT 



40 




Fig. 18. Error correcting performance comparison of the proposed QC LDPC codes, the PEG LDPC codes, and the PEG QC 
LDPC codes. 



VII. Conclusions 

The subgraphs of protographs, which cause inevitable cycles in the QC LDPC codes, are fully 
investigated in allowance with multiple edges through the graph-theoretic approach. For regular 
QC LDPC codes with girth larger than or equal to 12, we propose a systematic construction 
method of protographs which avoid inevitable cycles of length less than 12 by using balanced 
ternary designs. For regular QC LDPC codes with girth larger than or equal to 14, we provide 
construction methods of all J x L protographs with column-weight three and the number 
of double edges J — 2 by using various block designs. These construction methods can be 
extended to construct regular protographs with smaller number of double edges and with column- 
weight larger than three. Also, a construction algorithm of QC LDPC codes from the proposed 
protographs is provided based on the work in [7]. To check the validity of the proposed QC 
LDPC codes, we show that the proposed QC LDPC codes have larger upper bounds on the 
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minimum Hamming distance than the QC LDPC codes Ufted from single-edge protographs. 
Finally, the error correcting performance of the proposed QC LDPC codes is compared with 
those of PEG LDPC codes and PEG QC LDPC codes via numerical analysis. 
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